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Certain classical statistical systems with strong local constraints are known to exhibit Coulomb 
phases, where long-range correlation functions have power-law forms. Continuous transitions from 
these into ordered phases cannot be described by a naiVe application of the Landau-Ginzburg- Wilson 
theory, since neither phase is thermally disordered. We present an alternative approach to a critical 
theory for such systems, based on a mapping to a quantum problem in one fewer spatial dimensions. 
We apply this method to spin ice, a magnetic material with geometrical frustration, which exhibits 
a Coulomb phase and a continuous transition to an ordered state in the presence of a magnetic field 
applied in the [100] direction. 
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I. INTRODUCTION 

Classical statistical systems with Coulomb phases, 
where long-range correlation functions have power-law 
forms, have been of great interest in recent years. Of par- 
ticular interest ar e con tinuous transitions from Coulomb 
to ordered phases, 'i^^is] -^Jiich are prima facie incompat- 
ible with the standard Landau-Ginzburg- Wilson theory, 
since neither phase is simply thermally disordered. 

Various classical systems exhibiting power-law corre- 
lation functi ons- — such as geometrically frustrated anti- 
ferromagnets,!^ crystalline water icd^Eiand close-packed 
dimers on square and cubic latticed — have in common 
local constraints arising from dominant terms in the mi- 
croscopic Hamiltonian. In an effective description of the 
Coulomb phase, the local degrees of freedom are replaced 
by a continuum field, with the constraints replaced by a 
requirement that the field be diver genceless. The power- 
law behaviour is then unders tood i n terms of the corre- 
lations of this solenoidal field.'SlSliSl 

This coarse-grained picture does not, however, allow a 
description of certain ordered phases where the discrete 
nature of the local degrees of freedom is crucial, and 
therefore cannot describe transitions into these phases. 
These include states with saturated magnetic moments, 
as well as ordered phases in dimer models with hard-core 
repulsion between dimers.I^ 

Here we will describe an alternative approach, which 
leads to an effective critical theory for continuous tran- 
sitions from Coulomb phases to such ordered phases. It 
is based on the standard mapping between a thermal 
phase transition in d spatial dimensions and a quantum 
phase transition at zero temperature in d — 1 dimensions. 
The local degrees of freedom are represented by hard- 
core bosons, allowing their discrete nature to be treated 
exactly, while the power-law correlations in the Coulomb 
phase are reproduced by the fluctuations of the Goldstone 
mode that appears when the bosons condense. 

One clear disadvantage of the method is that it nec- 
essarily treats the system anisotropically, with one of 
the spatial dimensions being mapped to the quantum 



imaginary time, but it should be noted that space-time 
isotropy is often restored in the continuum limit of vari- 
ous quantum models. 

In a forthcoming work, we will present the applica- 
tion of this approach to the ordering transition of close- 
packed dimers o n th e cubic lattice .l^ Here, we address a 
model of spin ic d^ ^ ^^ l in an applied magnetic field aligned 
along the [100] crystallographic direction.'^ This system 
is convenient firstly because the mapping to the quan- 
tum problem takes a particularly transparent form, and 
secondly because the isotropy of the classical system is 
already broken by the magnetic field. The problem is 
also interesting in its own right because the phase tran- 
sition is produced by an external field, and because, as 
we explain in Sectio n |IC| it is an unusual example of a 
Kasteleyn transition^SEHTin three dimensions. 

A brief account of the application of a simplified ver- 
sion of this method to spin ice has appeared elsewhere,'^' 
along with results from Monte Carlo simulations. It 
was shown there that a [100] magnetic field can cause 
a Kasteleyn transition from the Coulomb phase to an or- 
dered phase where the spins are aligned with the field 
and fluctuations are frozen. For the corresponding quan- 
tum problem, this becomes a quantum phase transition 
from a Bose condensate to a vacuum state. 

Here we present the method in full, treating properly 
the crystal structure of spin ice, and show that doing 
so leads to an effective quantum Hamiltonian with non- 
hermitian directed-hopping terms. We demonstrate that 
these terms are crucial to understanding the Coulomb 
phase correlation functions via the quantum mapping. 
At the transition, however, we show they have no effect, 
so that the standard critical theory for the vacuum tran- 
sition of two-dimensional bosons applies, as conjectured 
in Ref. d 

In the remainder of this section, we will present the 
Hamiltonian for spin ice and review how, in the absence 
of an applied magnetic field, the correlation functions in 
the Coulomb phase can be understood in terms of an ef- 
fective solenoidal field. In Section|ll] we describe in detail 
the mapping to a quantum problem and the general form 
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FIG. 1: A section of the pyrochlore lattice, with the [100] 
crystallographic direction shown vertically. The sites of the 
lattice are at the vertices of the tetrahedra and the two (100) 
planes shown are one possible choice to define the mapping to 
a quantum problem described in Sectionjll] The points shown 
on these planes indicate the lattice sites for the quantum 
bosons, with black and white points respectively indicating 
the two distinct sublattices. We define the [Oil] direction as 
parallel to the bold lines drawn diagonally across the planes, 
connecting nearest neighbour sites. 

of the Hamiltonian that results, while in Section [IIT| we 
find the continuum action that corresponds to this Hamil- 
tonian. This allows us to obtain the critical theory for 
the ordering transition and also to show that the quan- 
tum model correctly reproduces the correlation functions 
within the Coulomb phase. Finally, in Section |IV| we 
show that these same correlation functions can be found 
directly from the microscopic quantum Hamiltonian, by 
treating the hard-core bosons as S = ^ spins and taking 
the lowest order in a spin-wave expansion. 



A. Spin ice: Hamiltonian 

The spin ice compounds consist of ions with large 
magnetic moments ('^ 10 Bohr magnetons) arranged on 
the sites z of a pyrochlore lattice (see Figure [l]) , with a 
strong crystal field favouring alignment of the spins along 
the lines joining the centres of neighbouring tetrahedra.^ 
At low temperatures, each moment Si is therefore con- 
strained to take the values = S'^e,, where Si — ±1 
is an effective classical Ising spin and is a unit vector 
joining the centres of two neighbouring tetrahedra. (The 
centres of the tetrahedra form a diamond lattice, and 
we choose the convention where always points from a 
certain diamond sublattice to the other.) 

In the spin ice materials Ho2Ti207 and Dy2Ti207, the 



largest interaction terms are long-range dipolar forces be- 
tween the magnetic moments.l^It has been shown,'^ how- 
ever, that the low-energy states of dipolar spin ice have a 
'projective equivalence' to those of a model with purely 
short-range ferromagnetic interactions, and we will work 
with this simplified model. Despite being ferromagnetic, 
this model is frustrated, and the ground-state degener- 
acy, determined experimentally by low-temperature mea- 
surements of the entropy,'i2l is in fact identical to that in 
the cubic phase of water ice. 

We treat spin ice in the presence of a field h along the 
[100] crystallographic direction of the pyrochlore lattice 
(the vertical direction in Figure [ij, whose effect is much 
smaller than the interaction between the spins. The field 
then acts as a perturbation within the manifold of degen- 
erate ground states of the h — problem. The Hamilto- 
nian can be written 

£ = - J ^ S,; • Sj - /in[ioo] • ^ Sj 
= JcS ^ SiSj — hcs ^(— l)^*S'i , 

where denotes that the sum is over nearest- 

neighbour sites i and j, and n[ioo] is a unit vector in the 
[100] direction, which we take as the z-axis. In the second 
equation, we have made use of the constraint Si — SiBi, 
and Zi denotes the z-component of the position of site 
i. Note that, in terms of the Ising degrees of freedrom 
Si, the effective coupling J^ff = | J is antiferromagnetic 
and heff = is an effective staggered field. Through- 
out, we will consider low temperature T and weak field, 
taking h,T J. 

For h — the energy is minimized by any configura- 
tion where every tetrahedron of the pyrochlore lattice has 
two spins pointing inwards and two pointing outwards; 
such configurations are said to obey the 'ice rules'.'^ The 
number of such states grows exponentially with system 
volume, and so there is no ordering even for T <ti J. 
The long-distance spin-spin correlation functions have 
a power-law form, which can be found using a coarse- 
grained solenoidal field, as we will review in Section [I B[ 

In the limit of large h, the Ising spins take the values 
that maximize the projection of the magnetic moment 
along the [100] direction. This results in the 'top' two 
spins in every tetrahedron (where [100] is defined as up- 
ward) pointing outwards, while the 'bottom' two point 
inwards. This is consistent with the ice rules and hence 
this state is the exact ground state for h > 0. 

As argued in Ref. [3 there is an ordering transition for 
T ~ /i in the limit of large J, with the ice rules apply- 
ing on both sides. One cannot straightforwardly write 
down a critical action of the Landau-Ginzburg- Wilson 
type for this strongly constrained transition.'^ Instead, 
we use an approach that relies on a mapping from this 
classical statistical problem to a quantum problem, which 
we describe in Section ITTl 
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Away from the infinite J/T limit, the ice rules are bro- 
ken on a finite density of tetrahedra by thermal fluctua- 
tions. In this case, the phase transition is replaced by a 
sharp crossover, which becomes more rounded off as J/T 
is increased (see the numerical results of Ref. |2] and in 
particular Fig. 2). 



B. Coulomb phase 

We first review a method for describing the behaviour 
on the disordered side of the transition .El We consider a 
coarse-grained version of the theory, where the spins 
on lattice sites i are replaced by a continuum field B(r), 
defined for all points in space r. While this gives the ap- 
propriate behaviour deep within the disordered phase, it 
does not incorporate the fixed spin length and so cannot 
be used to describe the transition to the ordered state. 

The ice rules require that the sum of over sites in any 
tetrahedron should vanish, which leads in the continuum 
limit to the constraint that B should be divcrgcncclcss: 
V • B = 0. Defining the gauge field A(r) by B = V x A, 
the only allowed continuum action respecting both the 
spatial symmetries and gauge invariance is 



St 



d^r (V X A) 



(2) 



where the omitted terms (involving higher derivatives) 
are irrelevant in the renormalization group sense. 

The correlation functions of the field B are therefore 
given by the standard dipolar correlation functions. 



(i3.(k)i3,(k)) ^ 6,, 



in momentum space and hence 



(B,(r)B,(0))^^':^'lJ^ 



(3) 



(4) 



in position space. The long-range correlation functions 
of the spins Si have identical forms. 

To determine the correlations of the effective Ising 
spins Si, we must take account of the four-site basis of 
the pyrochlore lattice. The correlation function decreases 
with increasing separation r as |r|~'^ and depends in sign 
and magnitude on the orientation of r and the choice of 
sites from the basis. As we explain below, we focus on 
the subset of sites that lie on a particular set of (100) 
planes illustrated in Figure [l] The sites lying in such 
planes are shown as the points in the figure, and form 
two sublattices, which we label a (black) and b (white). 

The directions along which the orientations of these 
spins are constrained are specified by the unit vectors Ga 
and eb, given by 



cos 9 



Ga =n[ioo] smW-hn[oii] 
eb = n[ioo] sin 9 - fifonj cos 



(5) 
(6) 



where tan6' = l/-\/2 and n^on] is a unit vector in the 
[Oil] direction, which is chosen parallel to the line joining 
neighbouring a- and b-sites. The spin Sg at a site on 
sublattice a can be written as Sg = 8363 (and similarly 
for b). 

We therefore find that the long-range spin-spin cor- 
relation function (in momentum space) for two spins on 
sublattice a is given by 



{S3S3 



(^[100] 



[Oil] 



(7) 



while that for two spins on opposite sublattices is 



{S3S,) ^ 



(In both cases, we have omitted a constant term, which 
upon Fourier transforming to position space becomes a 
delta function at the origin, and therefore has no effect 
on long-distance properties.) 

Note that the correlation function for sites within the 
a-sublattice is maximal along a line A;[ioo] = V^k^oii] in 
momentum space, parallel to the line (in real space) join- 
ing the centres of tetrahedra that meet at a-sites. By 
contrast, the correlation function between a- and b-sites 
has no such directionality, being symmetric in ±fc[oii]. 

In Section |IIIB| we will show that the same form for 
the correlation functions results from using the mapping 
to a quantum problem: see Eqs. (20 1 to (22). 



C. String picture 

We now consider the ordered phase, where the coarse- 
grained description is no longer appropriate, and the 
transition from this phase into the Coulomb phase. For 
T h J , all moments are as aligned with the applied 
field as possible, given the strong crystal field. As noted 
above, this implies that the top two spins in every tetra- 
hedron point out and the bottom two point in, as shown 
on the left-hand side of Figure |2] As T is increased, ex- 
citations above this ground state will start to appear, 
leading to a transition into the Coulomb phase described 
in Section |IB[ For T,h <C J, however, a single fiipped 
spin is not a low-energy excitation: it breaks the ice-rule 
constraint, and hence costs an energy on the order of J. 

In order to satisfy the ice rules, every tetrahedron must 
have the same number of fiipped spins at the top as at 
the bottom; an example is shown on the right-hand side 
of Figure [2] This in turn implies that every (100) plane of 
the lattice should have the same number of fiipped spins, 
and the lowest-energy excitation consistent with the ice 
rules is a string of fiipped spins spanning the system in 
the [100] direction. Such a string has a Zeeman energy 
of 2h/y/3 per unit of its length. 

An isolated string can be modelled as a random walk 
through the lattice in the [100] direction, with two choices 
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In the infinite-J limit, the description in terms of 
strings remains vaUd for all h: any state obeying the 
ice rules can be constructed by starting from the h ^ oo 
configuration and flipping all spins along a set of strings, 
with each spanning the system in the [100] direction. In 
other words, there is a one-to-one mapping from spin con- 
figurations obeying the ice rules to string configurations P 



II. BOSON WORLD-LINES 



FIG. 2: Two possible arrangements of the spins on a single 
tetrahedron of spin ice, where, as described in Section [1 A| a 
strong crystal field constrains the spins to point either directly 
towards or away from the centre of each tetrahedron. Both 
configurations have two spins pointing outwards and two in- 
wards and so minimize the ferromagnetic interaction energy; 
they are said to obey the 'ice rules'. The arrangement on the 
left shows the lowest-energy configuration in the presence of 
an applied magnetic field in the [100] direction (chosen ver- 
tically upwards, as in Figure [T]), with the spins maximally 
aligned along the field. On the right-hand side, one of the 
first excited states is shown, with two spins flipped with re- 
spect to the ground state. In order to obey the ice rules, there 
must be the same number of flipped spins at the top of the 
tetrahedron as at the bottom, and so the number of flipped 
spins is the same in each (100) plane of the pyrochlore lattice. 



at every step. There is therefore an entropy of In 2 per 
unit length of the string, and the free energy F in the 
presence of a single string ifl 



V3 



h-Tlii2 



(9) 



where Lz is the length of the system in the z-direction. 
As a result, for — > oo, strings are absent below a 
critical temperature Tk — 2/i/\/31n2. For T > Tk, how- 
ever, the entropic gain from introducing these string ex- 
citations into the system outweighs their Zeeman energy 
cost, giving a phase transition into the Coulomb phase. 

As noted previously,^ this transition has a strongly 
asymmetr ic cha racter and is an example of a Kasteleyn 
transition.l-LS'i^On the ordered side of the transition, fluc- 
tuations are completely suppressed and the transition ap- 
pears to be of first order, while there are fluctuations in 
the Coulomb phase and a divergent correlation length as 
the transition approaches, as for a continuous transition. 

The entropic argument as presented so far is only valid 
in the limit of inflnite J, where the ice rules are strictly 
enforced on every tetrahedron. Finite J/T allows thermal 
excitations that break the ice rules at a small density of 
tetrahedra, leading to a characteristic string length ^ 
e"'^/"^, where a is a numerical constant. For Lc < L^, 
this characteristic length replaces Lz in Eq. ([9]). The 
effect of finite J/T is therefore similar to the effect of 
a finite system size, replacing the phase transition by a 
crossover, whose sharpness increases with J/T. 



In order to understand the phase transition in detail, it 
is convenient to treat the strings as world-lines for bosons 
moving in a two-dimensional space, with the [100] direc- 
tion taken as imaginary time. The classical partition 
function is then interpreted as the partition function for 
a quantum problem, with the length of the system in 
the imaginary-time direction corresponding to the inverse 
temperature. In the thermodynamic limit we therefore 
have a two-dimensional quantum problem at zero tem- 
perature. 

Since the strings span the system in the [100] direction, 
the boson number is conserved and there are no closed 
loops, which would correspond to particle-antiparticle 
pairs. The mapping therefore provides a way to view 
the Kasteleyn transition as the zero-temperature quan- 
tum phase transition between the vacuum and a Bose- 
Einstein condensate, as the chemical potential for bosons 
fj, increases through some critical value fiQ. 

The vacuum of bosons then corresponds to the satu- 
rated state of the magnet at high field or low temper- 
ature. The lowest excitations in the quantum problem 
are single particles, which correspond to isolated strings. 
These have finite energy above the ground state and so 
the correlation functions are not long ranged. In fact, 
since the classical observables map to operators that are 
diagonal in the boson-number basis, their connected cor- 
relation functions are strictly zero in this phase. 

In the other phase, with a Bose condensate, there is 
a Goldstone mode corresponding to the broken phase- 
rotation symmetry, leading to power-law correlation 
functions for the spins. We will show in Section |III| 
that the correlations have the expected dipolar (three- 
dimensional) spatial dependence (a simplified treatment 
has been outlined previouslj^) . These long-range corre- 
lations remain as the field is decreased to zero, when they 
agree with the predictions of Section |IB[ 



In Section II C we will provide further evidence for the 
identification of the Coulomb phase with the condensed 
phase of the bosons, by showing that the off-diagonal 
long-range order in this phase can be interpreted in terms 
of deconfinement in the Coulomb phase. 

As conjectured earlier,^ and as we will demonstrate 
explicitly in Section |III A| the condensation phase tran- 
sition has dynamical critical exponent z = 2, and so 
our (2 + l)-dimensional system is at its upper criti- 
cal dim ension. Thermodynamic results for the boson 
problerrP^IiSIlI! can then be carried over directly to the 
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spin systemP 

As an aside, we note that the effect of finite J/T can 
also be straightforwardly understood using this picture. 
Tetrahedra that break the ice rules, which appear with 
nonzero density away from the infinite- J limit, become 
the ends of strings. In the quantum Hamiltonian, these 
are generated by source terms for the boson operator, 
which break the phase rotation symmetry and eliminate 
the phase transition. As noted previously, a large but 
finite J/T therefore has the effect of replacing the phase 
transition with a sharp crossover. 
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A. Effective Hamiltonian 

We will now present a more detailed analysis of the 
mapping from the problem defined in terms of spins to a 
quantum problem of bosons, which will allow the calcu- 
lation of correlation functions, for both zero and nonzero 
field. The mapping essentially follows the standard pro- 
cedure of using a transfer matrix to connect the degrees 
of freedom in one layer of the system to those in the next, 
followed by interpretation of the transfer matrix as the 
exponential of a quantum Hamiltonian. However, care 
must be taken to account properly for the pyrochlore 
structure. 

We first divide the lattice sites by their z coordinates, 
which we take to be measured along the [100] direction. 
Each spin in the pyrochlore lattice has neighbours only 
in the same plane and in those immediately above and 
below. The partition function can therefore be written in 
terms of a transfer matrix whose rows and columns are 
labelled by the possible configurations of a given plane of 
the lattice. 

The structure of pyrochlore is such that the positions 
of the sites in neighbouring (100) planes do not coin- 
cide, and it is in fact necessary to translate by four lay- 
ers before the lattice structure in a plane repeats. This 
is illustrated in Figure [3] which shows the projection of 
the pyrochlore lattice onto the (Oil) plane. The thick 
horizontal lines show the positions of two nearest planes 
with the same site positions, which are separated by three 
planes with sites in different positions. We therefore de- 
fine a new transfer matrix corresponding to the product 
of four elementary transfer matrices, which relates two 
such planes with the same lattice structure. 

This transfer matrix T can then be interpreted in the 
standard way as the exponential of the Hamiltonian ma- 



trix 7i, by writing T = e 



-•HAt 



where At is a step in 



imaginary time. The trace over a product of transfer 
matrices corresponds to the standard sum over histo- 
ries giving the partition function for the two-dimensional 
quantum problem. 

While it is possible to construct the matrix T explicitly 
by considering all arrangements of the spins in the three 
intermediate layers, and therefore in principle possible to 
find Ti exactly, we do not expect this to be necessary for 
our purposes. Instead, we will attempt to understand the 




[011] 



FIG. 3: Projection of the pyrochlore lattice onto the (Oil) 
plane. The sites of the pyrochlore lattice are at the vertices. 
The tetrahedra appear as triangles, and the dashed lines indi- 
cate tetrahedra that lie out of the plane. The two thick hor- 
izontal lines correspond to the thick lines shown in Figure [l] 
and represent one choice for the (100) planes that define the 
quantum problem. The distance between them in the [100] 
direction is one step in imaginary time, Ar. The points show 
a possible path for a string between these two planes, and the 
arrows indicate all possible locations for the end of the string, 
given its starting position. There are more points to the left 
of the starting location than to the right, and so a boson at 
this point has a preference for leftward hopping. 



general form of the Hamiltonian based on the symmetries 
of the problem, while restricting to the simplest version 
that is necessary to describe the large-scale behaviour 
correctly. 

First note that, since it describes a classical statistical 
model, the matrix elements of T are required to be real 
and nonnegative. We can therefore restrict to the case 
where TL is also a real matrix. However, as is often the 
case for classical statistical problems, the transfer matrix, 
and hence Hamiltonian, are not required to be hermitian. 
In fact, we will find (see Sections IIIB and IV C I that 
the sublattice structure of the correlation functions (see 
Section I B I is crucially dependent on the nonhermitian 



terms in the Hamiltonian, and so we will pay particular 
attention to such terms. 

The fact that T has all positive elements implies that 
its eigenvalue with largest real part — corresponding to 
the ground-state energy of 7i — is unique (within each 
boson-number sector), is positive, and has an eigenvec- 
tor with all positive weights. This can be interpreted as 
meaning that the bosonic Hamiltonian TL is unfrustrated, 
which is consistent with our assumption of a superfluid 
ground state. 

As described above, we will represent the state of a 
given layer by comparing to a reference configuration in 
which all spins have a positive projection on the +z axis. 
This configuration will be referred to as the vacuum of 
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bosons, while flipped spins are mapped to the presence 
of a boson. The bosons therefore have on-site hard-core 
repulsion by construction, and as noted above, the ice 
rules enforce conservation of boson number. 

We define the bosonic annihilation and creation oper- 
ators bi and bl for sites i in the two-dimensional lattice. 
The hard-core constraint on the bosons is implemented 
by the operator identities bibi — blbl — 0, so that the 
number operator rij — b\bi is restricted to the values 
and 1. (Operators for different sites commute.) 

We will divide the Hamiltonian Ti. into potential terms 
V and kinetic terms K, by writing 

n^V + K, (10) 

and address these two parts in turn. 

The simplest contribution to F is a chemical potential 
term Vq giving the energy cost of adding a boson world- 
line to the system. In the presence of an applied magnetic 
field along the +z axis, every flipped spin has Zeeman 
energy cost 2/icff, so that we can write 

where ^ oc —h. The classical model is symmetric under 
h ^ —h along with an inversion of all the spins. After 
mapping to the quantum problem, this takes /i —fi 
and exchanges particles and holes; the constant (— |) has 
been chosen to make this symmetry explicit. A vanishing 
applied field corresponds to half filling in the boson rep- 
resentation, and at this point the model is particle-hole 
symmetric. 

Interactions between the bosons are dominated by the 
hard-core repulsion, but there will also be off-site inter- 
actions that contribute to V. These terms are necessarily 
hermitian and we will not consider them further. 

The kinetic terms K will include both quadratic hop- 
ping and more complicated correlated hopping terms, 
and, as we have observed above, are not required to be 
hermitian. For a single boson, the physical origin of the 
nonhermitian terms is clear from the lattice structure. As 
shown in Figure [3] a string starting on a given site has 
a tendency to hop in one direction in preference to the 
opposite direction, and it is this directed hopping that 
requires the Hamiltonian to be nonhermitian. 

Note, however, that a hole surrounded by particles 
hops preferentially in the same direction as an isolated 
particle, as can be checked by exchanging particles and 
holes in Figure [3] As a result, the quadratic hopping 
term, 

K, = -Y,t^AbJ, (12) 

is required to be symmetric under b^ so that tij — 

tji and Ki is hermitian. The same conclusion follows 
from the requirement of particle-hole symmetry when 
/i = 0. 



We must therefore include higher-order terms in order 
to describe the directed hopping, and so we seek cor- 
related hopping terms that are particle-hole symmetric. 
The simplest such term is 

K2 = -Y,w,,,{nt-\)b\b^, (13) 

ijt 

where Wiji = —Wja and so K2 — —K\. Note that near 
the transition, where the density is small, K2 reduces to 
a quadratic hopping term of the form of Ki, but with 
effective directed hopping = — 1 '^^Wijt. 

We therefore take V = snd K = Ki + K2. The 
coefficients tij and Wiji are required to be real and chosen 
so that, as noted above, the hopping is unfrustrated. 

B. Symmetries 

In constraining the general form of terms that can ap- 
pear in the effective Hamiltonian and the corresponding 
continuum action, it is important to consider symme- 
tries inherited from the original problem defined on the 
pyrochlore lattice. The choice of z as the imaginary time 
axis and of the particular (100) planes as described in 
Section [ll A| reduces the symmetry group somewhat. 

This is illustrated in Figure |4] which shows the py- 
rochlore lattice projected onto the (100) plane. It can 
be seen immediately from the arrangement of tetrahedra 
that nearest neighbours in the x-direction are equivalent, 
whereas those in the y-direction are not. The lattice can 
therefore be represented as two interpenetrating square 
lattices, indicated by the filled and empty circles. We 
define the elementary translation operators ICx and ICy 
which take sites into their nearest equivalent neighbours 
in the x- and y-directions respectively, and which are 
symmetries of the Hamiltonian. (Note that K-x and ICy in 
fact translate by the same absolute distance, although JCx 
connects neighbours, while ICy connects next-neighbours 
in the y-dircction.) We also define the operation JCz, a 
translation by four layers in the z-dircction, which cor- 
responds to time-translation symmetry in the quantum 
problem. 

One can in addition see that there is a reflection sym- 
metry, which we label TZx, about the red line labelled a, 
which takes x ~x but does not exchange sublattices. 
There is similarly a reflection symmetry TZy about the 
green line labelled b, taking y — > —y, but in this case the 
symmetry exchanges the two sublattices. 

Besides TZx and TZy, there is a third reflection symme- 
try, which involves reflection in the blue line labelled c, 
and also takes y — > — y. This does not affect the sublat- 
tices, but changes the height of the tetrahedra: the tetra- 
hedra drawn with thick lines in the figure are exchanged 
with those drawn with dotted lines and those with thin 
lines are exchanged with those with dashed lines. This 
amounts to a reflection in the x-y plane (i.e., z —z), 
which taken together with the reflection about the x- 
direction constitutes a third reflection symmetry, Xy. In 
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FIG. 4: Projection of the pyrochlore lattice onto the (100) 
plane, with the x- and y-axes corresponding to the [Oil] and 
[Oil] directions respectively. Each tetrahedron is shown as 
a square with a line joining the two uppermost points. The 
weight of the lines indicates the height of a given tetrahedron 
in the [100] direction, with thick lines highest, followed by 
thin, dashed, and dotted lines respectively. The strucure is 
periodic, so that another tetrahedron of the same orientation 
lies four layers below (and above) each tetrahedron that is 
shown. The circles indicate sites in a particular (100) plane, 
which, after the mapping to the quantum problem, become 
the sites for the bosons. Filled and empty circles correspond 
to the black and white points shown in Figure [T] and denote 
the two sublattices, which are distinguished by the structure 
of the surrounding lattice. The thick red and green lines (la- 
belled a and b) show planes of reflection symmetry, while 
there is a symmetry under rotation by 180° about the blue 
line (labelled c). Rotation by ^ about the cross (labelled d) 
and translation in the [100] direction maps sites in one (100) 
plane into those in the plane immediately above. 



three dimensions, ly is a rotation by 180° about the x- 
direction, while, after mapping to the quantum problem 
defined in two dimensions, it becomes a ^/-reflection ac- 
companied by a time-reversal operation. 

Finally, consider once more the elementary transfer 
matrix linking one (100) plane to that immediately above 
it. As explained above, this matrix does not act within 
the quantum Hilbert space, since it relates sites on two 
different (two-dimensional) lattice structures. Instead, 
we define the operator given by this elementary 

'evolution' operator followed by a rotation by | about 
the cross labelled d in Figure |4] This composite trans- 
formation returns the sites to their original locations and 
hence acts within the Hilbert space. It is also clear that it 



commutes with the Hamiltonian, since the original choice 
of a particular (100) plane was arbitrary. 

As for the full transfer matrix T, we will not write 
down the precise form of the operator 7^/4, and we will 
not attempt to find the general constraints it implies for 
the Hamiltonian. Instead, we will merely note that, in 
a coarse-grained description, it has a particularly sim- 
ple effect on operators that do not distinguish the two 
sublattices, giving only a rotation in the x-y plane. 

C. Off-diagonal long-range order and 
deconfinement 

As we have explained above, we identify the Coulomb 
phase of the classical statistical model with the condensed 
phase of the quantum bosons. In Section |IIIB| we will 
show that the long-range correlations in the quantum 
model reproduce those of the classical model, support- 
ing this identification. Here, we present a qualitative ar- 
gument based on the concept of off-diagonal long-range 
order (ODLRO). 

First, note that while the quantum model has {bi) 7^ 
in the condensed phase if total boson number is allowed 
to fluctuate, there is no equivalent order parameter in the 
Coulomb phase. Under the mapping we use, boson num- 
ber is strictly conserved and a nonzero order parameter 
{hi) is impossible. Spontaneous symmetry breaking in the 
quantum problem is therefore hidden in the statistical- 
mechanics equivalent. We instead consider ODLRO: in 
the condensed phase, the correlation function {b\h^) has 
a nonzero limit for large separation of the sites i and j, 
given by the square modulus of the order parameter. 

Translating this quantum expectation value into the 
language of the statistical problem gives {b\bj) ^ Zij/ Z. 
Here, with (3 the quantum inverse temperature (to be 
taken to infinity), Z = TrT'' is the partition function, 
and Zij is given by 

Z,,^TrTH\b^. (14) 

This can be interpreted as the partition function, calcu- 
lated in the presence of one string (boson world-line) that 
ends at the point j on a certain arbitrary (100) plane and 
another that begins at the point i on the same plane. 

This is most simply interpreted by returning to the 
representation in terms of solenoidal fields described in 
Section |IB[ A string ending at a given point corresponds 
to a tetrahedron on which the ice rules are broken, so 
that the coarse-grained field B(r) has nonzero divergence. 
We therefore have two monopoles of opposite 'charge' 
separated by a large distance. 

In a confining phase of the field B, the energy in- 
creases linearly with separation, so that Zij/Z is ex- 
ponentially suppressed for large separations. In the 
Coulomb phase, however, monopoles are deconfined, and 
Zij/Z approaches a constant, given by the monopole fu- 
gacity, for large separation. 
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The ODLRO in the quantum superfluid is therefore 
equivalent to monopole deconfinementfiS in agreement 
with our identification of the Coulomb and superfluid 
phases. 



III. CONTINUUM THEORY 

In this section, we show that the long-range form of 
the correlation functions in spin ice, with and without 
an applied field, can be found using an effective con- 
tinuum theory for the hard-core bosons. We also find 
a critical theory for the Kasteleyn transition, and show 
that, although the directed hopping is necessary for the 
appropriate long-range correlation functions, it does not 
affect the critical behavior. In Section |IV| we show that 
the same long-range behavior can be found directly from 
the microscopic quantum Hamiltonian, by using a large- 
S expansion. 

Our approach here is to start with the Hamiltonian 
Ti. and write down the most general continuum action 
consistent with the symmetries of the hard-core boson 
problem, described in Section |IIB[ In the Appendix, we 
sketch an approach that would in principle allow such a 
continuum action to be derived directly. The action is 
expressed in terms of bosonic fields corresponding to the 
hard-core boson operator b, in the limit where the spatial 
coordinates x and y and the imaginary time t are taken 
as continuous. Since it is important to preserve the two- 
sublattice structure, we define (c-number) fields tpa and 
-013 corresponding to the two sublattices. 

Besides the spatial symmetries, there is a phase- 
rotation symmetry arising from the conservation of the 



boson number, which leaves the action invariant un- 
der the simultaneous transformations ipa —> V'ae'^ and 
V'b — *■ V'bG'^- Note that only the total number of bosons 
on both sublattices is conserved, and so terms such as 
0*'0[j are allowed by this symmetry. 

We now consider each of the symmetries defined in 
Section II B in turn. The translations ICx, ICy, and ICz 
shift X, y, and r (infinitesimally, in the continuum limit) 
and so terms with explicit dependence on the coordinates 
are forbidden. 

As for the reflection symmetries, TZ^ simply takes x — > 
—X, while TZy takes y — y, simultaneously exchanging 
the two sublattices. The derivative operator dx = d/dx 
can therefore only appear in the combination d^, but a 
single j/-derivative can appear in a combination such as 

In the continuum theory for standard lattice boson 
models one also generically has symmetry under time 
reversal accompanied by particle-hole exchange: dr 
—dr, tp —> tp*] and, at points where the theory is particle- 
hole symmetric, such as half-filling for hard-core bosons, 
the two transformations are separately symmetries. In 
our problem, the nonhermitian nature of the Hamiltonian 
means that this is no longer the case. Instead, as we show 
in the Appendix, the refiection operation Xy leads to a 

-dr, Va.b ^ Cb- At 



symmetry under dy —dy, dr 
half-filling, the particle-hole transformation "^a b ^ V'a b 
is again a separate symmetry, which we label Xph. 

The most general Lagrangian density that is consistent 
with these symmetry constraints, up to quadratic order 
in the fields and second derivatives (and omitting terms 
that are total derivatives), can be written 



+ \Ky{r.dyA - AdyA) + [{d^AmA) + {d^.A){dM + a; Kd^AKdM + {d^A){dM] 

I 



(15) 



where a sum over G {x, y, r} is implied. At half-filling, 
particle-hole symmetry implies Kr = k'^ = Ky = 0. (Note 
that we so far have made no use of the symmetry 7^/4, 
which has a nontrivial effect on the sublattice labels a 
and b.) 

To evaluate the correlation functions, it is also neces- 
sary to find expressions for the boson number operators 
in terms of the fields 4'a and Ai and symmetry consid- 
erations can again be used. Defining fields and nb 
representing the number fiuctuations on the two sublat- 



tices, we have nt, unuer /c^ anu na,b ~^ ~'"-3,b 

under Xph; they are unaffected by all other symmetries. 
(Note that these are the number fluctuations; the ac- 
tual number operators also contain constant terms, which 
will cancel when calculating connected correlation func- 
tions.) These symmetries are not, however, sufficient to 



under TZ„ and 



determine the number operators completely. We find, 
for instance, that Us + rib is given by a linear combina- 
tion of KAdrA + AdrA)' KAdrA + A9rA), and 
KA^v''Pa ~ '0b^i/V'b)i where the relative coefficients can- 
not be determined by this method. 



A. Critical theory 



Using the Lagrangian density given in Eq. (15 



can address the behaviour at the Kasteleyn transition 
from the vacuum to a phase with a Bose-Einstein con- 
densate. This necessarily occurs away from half filling, 
where there is no particle-hole symmetry (although an 
identical transition, related by a particle-hole transfor- 
mation, occurs at unit filling). At the transition, both 
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ips and tpt, acquire nonzero expectation values, breaking 
phase-rotation symmetry. 

It is convenient to transform from i/'a and ip^, to a new 
basis that diagonahzes the constant, quadratic part of 
C Due to the symmetry TZy, which exchanges the two 
sublattices, the appropriate fields are — ^paiipt,, and 
the constant part of the Lagrangian density then becomes 



(16) 



The difference of the two constants r+ and r_ is given by 
r2 in Eq. ( |15| , and will generically be nonzero. Both 
and r_ will be positive in the vacuum phase and will grow 
smaller as the transition is approached. At the transition, 
one of the two (strictly, its renormalized value) changes 
sign, and the corresponding field V'± acquires a nonzero 
expectation value. As a result, (tpa) and (ipb) both be- 
come nonzero, and with the same magnitude. 

The theory can then be rewritten in terms of the crit- 
ical field, V'+ say, while the other (■0-), which still has a 
'mass' term is integrated out. To determine the form of 
this critical theory, consider the quadratic terms with 
single derivatives. Rewriting the term in C with one 
dy derivative in terms of ^p± gives ipldyip^ — ip'^dy'ip^ ~ 
4'+dyip_ + 4'*-dy4'-\.- Integrating out the field ■0- leaves 
only a term of the form {dyip*j^){dy'ip_^). On the other 
hand (as long as Kr + 0); there will be a term 

'4'^drip_^_ in the Lagrangian density, which remains after 
ip- is integrated out. 

Furthermore, since the field "0+ has no dependence on 
the sublattice structure, we make use of the symmetry 
Ti/4 to require invariance under a ^ rotation in the x- 
y plane. (The critical theory therefore has an emergent 
symmetry under this operation, and in fact under con- 
tinuous rotations in the x-y plane.) The critical theory, 
which is written only in terms of , then takes the form 

-Ccritical = »'+V-'+'0++i0+9r0'++A+(9j0+)(aj-0+)+ ••■ , 

(17) 

where a sum over j € {x, y} is implied: the second time 
derivative has been dropped. This therefore describes the 
standard vacuum transition of bosons, with dynamical 
critical exponent z = 2, as conjectured previously.^ 



B. Condensed phase 



We can also use the continuum theory of Eq. ( 15 ) 



to calculate correlation functions within the condensed 
phase. Once the fields ips and ipt, acquire nonzero ex- 
pectation values, they are most conveniently described 
in terms of phase and amplitude modes, by writing 

0a = Vpo + SpJ^*+'^ and = v/A^+^e'(*-«) , 

(18) 

where po is the (real) condensate density, and 5pa,b, 0, 
and 9 are real fields. 

The amplitude modes, Sps and Spt,, and the relative 
phase mode, 9, correspond to gapped excitations, while 



the overall phase is a Goldstone mode, resulting from 
the broken phase-rotation symmetry of the condensate. 
Only this gapless mode can contribute to the long-range 
correlation functions, and we therefore integrate out the 
remaining gapped modes. While this process can be car- 
ried out explicitly, it is clear from symmetry considera- 
tions that the general form for the resulting Lagrangian 
is simply /Z^ = {dr(l)y + {djCp)^ + ■ ■ ■ , where a sum over 
j e {x, y} is again implied. The constant c, giving the 
speed of sound in the condensate, cannot be determined 
by symmetry. 

In terms of the field (p, the number operators take a 
simple form: symmetry requires that we have 

ria ^ {dr + vdy)(l) and rib ^ {dr ~ vdy)(p , (19) 

but does not fix the magnitude or sign of v. 

Taking the Fourier transform therefore gives for the 
correlation functions 



{n^nb) 



{bj + vkyY 



{lo — vky)^ 
iJ^ + c?(kl + fc2) 

ti)-^ — v^k^ 
io^ + c^{kl + kl) 



(20) 
(21) 
(22) 



After the identification of w with fc[ioo] and ky with fc^oii]; 
these become exactly equivalent to Eqs. ([t]) and ([s]), de- 
rived in Section |I B| using the mapping to a solenoidal 
field. Note that Eqs. ^ and ^ apply only in the ab- 
sence of an applied magnetic field, whereas Eqs. (20) to 



( |22| are valid throughout the Coulomb phase, with the 
dependence on h appearing in the constants c and v. 
These cannot be determined by this method, however, 
since the quantum mapping explicitly reduces the sym- 
metry of the three-dimensional problem by picking out 
the z-direction. 

As the vacuum transition is approached from the con- 
densed phase, the speed of sound decreases continu- 
ously towards zero, and standard results for the boson 
problenF^ give c '-^ p — pq. This leads in spin ice to 
anisotropy in the correlation functions between the z- 
direction and the x- and y-directions, with relative scale 
c ~ -y/ZiK — h, where /ik is the critical field at the Kaste- 
leyn transition. 



IV. LARGE-5 CALCULATION OF 
CORRELATION FUNCTIONS 

We now present an alternative calculation of the spin- 
spin correlation function in the Coulomb phase, starting 
from the Hamiltonian in the bosonic representation. Our 
approach is to rewrite the hard-core bosons in terms of 
5=1 spins and then use a spin-wave approximation 
based on a large- S* expansion. 
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A. Holstein-Primakoff expansion 

We begin with the effective Hamiltonian described in 
Section [tfll 

n = Vo + Ki+K2 

--^^i:{n.-'^-^UAb, (23) 

iji 

where the coefficients /i, tij and Wije, are real and obey 
tij — tji and Wij£ — —Wjie- We first map the hard- 
core bosons into 5 = spins by identifying Ui = blb^ = 

i-Sf , b] = and h = S+, where Ef = Ef ±iSf . (We 
use the symbol S rather than S to avoid confusion with 
the classical spins of the original spin ice Hamiltonian, 
to which these quantum spins are not simply related. In 
particular, the labels x, y, and z do not correspond to the 
directions of the classical spins or the crystal structure.) 
In terms of spins, the chemical potential becomes an ap- 
plied field along the z direction, while the kinetic terms 
produce couplings between the x and y components of 
nearby spins. 

We now extend this model to 5 > 5 , so that become 
spin-S' operators for general S. For large S, the Holstein- 
Primakoff transformation can be used to expand around a 
classical ground state. Such states have all spins parallel 
and at an angle 6 to the field, and we choose the state 
where they lie in the x-z plane. Rotating to new axes 
aligned with the classical spins, we write 

E^ = cos6'- S^sine* (24) 
^ ^^sm9 + t^cos0 (25) 
Sr = Sf. (26) 

The Holstein-Primakoff transformation is given by 




(27) 

(28) 
(29) 



where = Ef ±iEj^, and a.^ and a] are boson operators 
obeying [aj,aj] = Sij. (Note that these bosons are not 
simply related to the original hard-core bosons of the 
quantum problem.) 

For large S, the square roots can be expanded as power 
series in S~^. While such an expansion should clearly 
not be expected to give quantitatively accurate results 
for the case S = | , we will show that it already captures 
the expected physics at quadratic order in the operators 

and a] , even at half filling, where the hard-core nature 
of the bosons is most important. 



The classical ground state can be found by applying 
Eqs. (24 1 to (29 1 to the Hamiltonian TC and keeping only 



those terms in the expansion that contain no boson op- 
erators. This gives 



Ho = Nij,S cos 9- S^Y^ U-j SIT? 



- S^^w^je cose siii^ , (30) 

iji 

where N is the number of sites in the lattice. The final 
term vanishes because Wiji = —Wju. Let to = J^j^ij^ 
which can be shown to be independent of i using a com- 
bination of the transformations ICx , ICy and TZy , defined 
in Section |IIB| The minimum of Tio is given by choosing 



COS! 



^J■ 



2Sto 



(31) 



provided that ^2Sto < /i < 2StQ, or equivalently, that 
the chemical potential lies within the dispersion band. 
At the transition into the vacuum, the chemical potential 
goes down through the bottom of the band; within the 
vacuum phase, /i < —2StQ, and the minimum is instead 
given by cos 9—1. 

With this choice for 9, the terms in the expansion that 
are linear in boson operators cancel. The next terms are 
quadratic in the boson operators and, as usual for the 
Holstein-Primakoff expansion, include terms of the form 
aiGj and a,la.]j as well as number-conserving terms such 

as aja^. After some algebra, these quadratic terms can 
be written as 



«-5EW iktZ 



where 

Aij = 2StoSij - Stij{l + cos^ ( 



(32) 



(33) 



Bij — 2S^ cos^ 9 Wiji + sin^ ^ ^^(w^j — wiji) 

(34) 
(35) 
(36) 



Cij — St-ij sin 



Dij = -S"^ sin^ 6* ^(w£ij -I- wiji 



Note that Aij, Cij and Dij are symmetric in their indices, 
while Bij is antisymmetric. The nonhermitian nature 
of TL is therefore reflected in 7^2: taking the hermitian 
conjugate of 0.2 gives Bij 



-Bij and Dij 



B. Bogoliubov transformation 

The quadratic Hamiltonian 7^2 can be diagonalized us- 
ing the standard Bogoliubov transformation of bosons, in 
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which some care must be taken because 7^2 is nonhermi- 
tian. First define the column vector a so that ai = ai 
for 1 < i < iV and a, = a]_^ iov N + I < i < 2N . 



The commutation relations become [aj,a^] — rjij, where 
is a 2N x 27V diagonal matrix with +1 for the first 
TV elements on the diagonal, and —1 for the remainder. 
We transform from a and its conjugate cx^ to a and a. 
(which are not hermitian conjugates), by 



a = Ya and = aW . 

The matrices V and W 
but are instead related 
[cxi, cxj] — rjij, so that W = r]V~^r]. 



(37) 



are not hermitian conjugates, 
to preserve the commutator: 



The quadratic Hamiltonian 7^2 can be written as 



1 



(38) 



(where a matrix product is implied), so we must choose 
V such that H = V^^rjHV is a diagonal matrix, giving 



7^2 = -^a.riH.a . 



(39) 



Note that, since neither H nor jyH is hermitian, the ele- 
ments of H are not necessarily real. 

Nonetheless, it can be showrP^l that for every eigen- 
value ei of T7H, there is a corresponding eigenvalue — e^. 
After appropriately ordering the eigenvalues and eigen- 
vectors, one can therefore identify the elements on as an- 
nihilation {\ < i < N) and creation {N + \ < i < 2N) 
operators and write 



N 



H2 = E 



e,a,;a,; 



(40) 



where Ci are the eigenvalues with positive real parts (and 
a real constant has been dropped). Note that di is not 
the hermitian conjugate of di, and so that number op- 
erator didi is not hermitian. The commutation rela- 
tions nonetheless ensure that di and di act in a bona 
fide Fock space, and that the number operator has real, 
nonnegative-integer eigenvalues. 

In our case, where the Hamiltonian and hence 
are invariant under the translation operators IC^ and JCy, 
the eigenvectors can be assigned to momenta within the 
reduced Brillouin zone. This contains -^N points, and 
we therefore expect two bands, of which only the lower 
band will be important for the long-range properties with 
which we are concerned. Let Ok and d^ be the operators 
for a state in the lower band with (crystal) momentum 
k, and let ek be the corresponding energy. These lower- 
band excitations are the phonon modes of the condensate 
and will therefore have a linear spectrum. 



Ek — yc^k'^ + CyA;^ , 
where Cx and Cy are the phonon velocities. 



(41) 



In terms of the coefficients defined in Eqs. ( 33 1 to ( 36 1 
Cx and Cy are given by 



cl = 2Stosm^eY,xliC^o-AM) 

i 

cl = 2Sto sin^ J2 y^C,^ - 



(42) 



D 



jO) 



(43) 



While the full spectrum is not necessarily real, the as- 
sumption that the ground state of the system is a stable 
superfluid implies that the phonon velocities are real. 



C. Correlation functions 

The connected part of the spin-spin correlation func- 
tion in spin ice, Cc(r,r';z), for spins at points r and r' 
in their respective (100) planes and separated by z in 
the [100] direction, translates in the hard-core boson lan- 
guage to the time-ordered number-number correlation 
function. At inverse temperature /3 (which we will later 
take to infinity), the expectation value of the number 
operator at site i is 



= ^ Tr (r^7 



(44) 



where Z = TrT'^ is the partition function. Note that 
the operators here refer to the hard-core bosons 6;, rather 
than the Holstein-Primakoff bosons . The pairwise cor- 
relation function (including disconnected parts) is 



j i Tr (r'3--n,T-7i 



for z > 0, 



Tl^ln,) for z < 0. 



(45) 



(Note that the relative units for distance in the [100] 
direction and within the (100) plane are arbitrary, and 
so we have chosen to take the unit of z as four layers in 
the [100] direction.) 

First observe that exchanging the two sites i and j is, 
as expected, equivalent to taking z — > — z; we therefore 
focus on z > 0. 

We rewrite the correlation function using a spectral 
representation, by diagonalizing the transfer operator as 



T = J2\E)e-''{E\, 



(46) 



where E labels a complete set of eigenstates of the Hamil- 
tonian. This gives 



= i^(i?|n,|i?)£ 



-f3E 



(47) 



12 



and 



C{t,j;z) = ^ E {EK\E'){E'\n,\E)e-^^P-^^e-'''^ 



E.E' 



(48) 



which reduce, in the zero-temperature hmit, to 

(n,) - {EQ\n,\Eo) (49) 

and 

C(z, j; z) = ^ (i?o|n,: |£;o)e-(^-^'')^ , (50) 

E 

where E^ labels the ground state. The term in the sum 
with E ^ Eq is independent of z and cancels when the 
connected correlation function Cc is calculated: 

Cc{iJ;z) ^C{i,j;z) ~ {n,){nj) (51) 
= ^ {Eo\n,\E){E\n,\Eo)e-^''~''"^^ . (52) 

E^Eo 

The translation symmetry of the Hamiltonian means 
that we can choose a set of eigenstates \E) which are also 
eigenstates of the translation operators ICx and ICy. In 
particular, the eigenvalue of the translation operators in 
the ground state will be +1. We let ria and Uh denote the 
number operators for the sites belonging to each of the 
two sublattices in the unit cell at r = 0. The translation 
operators can then be used to give 



C,{i,j;z)^ J2 {EoME){E\n^^\Eo)e 

E^Eo 



-{E-E„)z^iKE 



(53) 

where € {a,b} denotes the sublattice that site i be- 
longs to, Tij is the separation of sites i and j within the 
(100) plane, z is their separation in the [100] direction, 
and Kb is the total momentum of the state E. 

The matrix elements {EqItIctIE) and {E\n„\EQ) (which 
are not complex conjugates, because of the nonhermitian 
Hamiltonian) are to be calculated consistently within 
the Holstein-Primakoff expansion. Using the mapping 
to quantum spins and Eqs. (24 1 to (29 1, we have 



Scose+^ + J^{a^+al)sme + 



(54) 



where the omitted terms are at least quadratic in the 
boson operators. The matrices V and W can then be 
used to express Ui and a] in terms of cik and Ok, using 



Eq. (37) 



The quadratic approximation to the Hamiltonian, TI2, 
has eigenstates with definite numbers of quasiparticle ex- 
citations, and the operator rii, which is linear in Ok and 
Sk, can only create or annihilate a single quasiparticle. 
The constant terms in the expansion for rii will there- 
fore always cancel when calculating off-diagonal matrix 
elements: these will involve only states E with a single 



quasiparticle. If the momentum of this quasiparticle is 
k, then K^; = k, and one can write 

{Eo\ia, + ai)\E) = /.(k) and {E\{a, + ai)\Eo) - /^(k) . 

The connected correlation function can then be written 
in the form of an integral over k. 



d^k 

(2^ 



/^rk)/' (k)e-^''^e 



(56) 

Note that, if and were complex conjugates, as 
would be the case for a hermitian Hamiltonian, this ex- 
pression would be symmetric under i *^ j. The nonher- 
mitian nature of the Hamiltonian is therefore crucial to 
the directional dependence of the correlation functions, 
as noted in Section III Al 

We have so far restricted to z > 0; the form for z < is 
given by exchanging i and j , or equivalently by exchang- 
ing Gi and (Tj and taking k — > — k. The Fourier transform 
to momentum- and frequency-space is therefore given by 



CT'(k,w) = ^sin 



/,(k)/;,(k) , 4,(-k)/;(-k) 



ek 



(57) 

As in the square-lattice case,'2lone finds /o-(k) ^ ^/jkl 
for small |k|. The anisotropy in the Hamiltonian leads, 
in this case, to a dependence on the direction of k of the 
form 



/a(k) 



(58) 



(and the same for /'), where ek is the dispersion given 
in Eq. (41 1, v is a real constant, and the sign of the 
imaginary part is -t- when a is in the a-sublattice and — 
for the b-sublattice. 

The characteristic velocity v can be written as 



V = 2Sto sin 



j:,{~l)v^{Aio + Cio) 



(59) 



in terms of the coefllcients defined in Eqs. (33 1 to (36) 



One therefore finds, for the connected correlation func- 
tion within the same sublattice. 



Caa(k,w) -sm e 2-^ 

2 U!^ 



2vkyU! 



(60) 



with V — > —V for Ct,b- For opposite sublattices we have 
instead 



S 

Cab(k,tj) ^ ^sin^ 



(61) 



It can immediately be seen that the former expression has 
a spatial asymmetry and hence a preferential direction, 
while the latter is symmetric in ±ky. 

These expressions are in fact exactly equivalent to 
those given in Eqs. ([7| and ^ in Section IB and in 
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Eqs. (20 1 to (22 1 in Section IIIB as can be shown by 
subtracting {Sj2) siv? 6 from both. (This term, which 
is independent of k and uj, becomes a delta function in 
real space and hence has no effect on the long-range be- 
haviour.) 



D. Intermediate asymptotics near transition 

As the phase transition is approached from the super- 
fluid, the chemical potential decreases towards the bot- 
tom of the band, and Eq. (31 1 implies cos 9 = 1 at the 



transition point. It follows that the phonon velocities Cx 
and Cy vanish like the square root of the chemical po- 
tential as the transition is approached, as noted in Sec- 



tion IIIB There is therefore an intermediate regime of 
behaviour within the superfluid but close to the transi- 
tion, where smO <C |k| ^ 1 and the linearized dispersion 



relation in Eq. (41 1 is no longer valid 



From Eqs. (|33|^to (|36|), one finds that Cij and Dij 



vanish as sin0 — > 0, and so the only terms that remain 
in the quadratic Hamiltonian TL2 are those that conserve 
(Holstein-Primakoff ) boson number. The Hamiltonian is 
therefore diagonalized with a straightforward canonical 
transformation, giving a dispersion relation of the form 



2m„ 



2m„ 



(62) 



phase at strong field was mapped to a quantum phase 
transition between a condensed phase of bosons and the 
vacuum. While one cannot naively write a Landau- 
Ginzburg- Wilson (LGW) theory for the original tran- 
sition, since neither state is thermally disordered, the 
quantum transition is described by th e standard critical 
theory for bosons at low density.tiHlS] 

The quantum Hamiltonian that results from the map- 
ping is strongly interacting, and we have not attempted 
to find it explicitly. Instead, we showed how symme- 
try considerations can be used to constrain the general 
form for the Hamiltonian and for the corresponding con- 
tinuum action. Using these general forms, we presented 
two different calculations of the correlations within the 
superfluid phase, which, in the zero-field limit, agree with 
previous results for the classical Coulomb phase.® ^ 

While the mapping from quantum mechanical systems 
in d — 1 dimensions to classical systems in d dimensions 
is of course standard, we believe that this is the first ex- 
ample of the use of the reverse mapping to derive a LGW 
theory for a naively non-LGW thermal phase transition. 
(See Ref. Jljfor an example of an application to an LGW 
transition in two spatial dimensions.) As noted above, a 
forthcoming work will present the application of this ap- 
proach to the ordering transition of close-packed dimers 
on the cubic lattice,'^ which is also believed to exhibit a 
non-LGW thermal phase transition. 



and constant matrix elements f^, f'^ ^ I + 0{\\i\). 

The connected correlation function is therefore given 

by 



ciKuj) 



(63) 



with no sublattice dependence, to leading order. After 
Fourier transforming back to real space, this gives 



C{r,z) - Aexp 



2H 



(64) 



the correlation function for a random walk in two 
dimensions .fSl This is to be expected near the transition 
into the vacuum, where the density is small and the corre- 
lations at distances shorter than the string separation are 
controlled by the behaviour of an isolated string. In this 
random-walk regime, the correlation function no longer 
refiects the directed hopping of the original Hamiltonian, 
since a single step is sufficient to randomize the sublattice 
of the string. 



V. DISCUSSION 

In summary, we have studied spin ice in a [100] mag- 
netic field by mapping the classical statistical system in 
three spatial dimensions into a quantum model in two 
spatial dimensions. The thermal phase transition be- 
tween a Coulomb phase at weak field and an ordered 
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APPENDIX: SYMMETRIES OF THE 
CONTINUUM THEORY 

Here we present a method by which the continuum 



action presented in Section HI could in principle be de- 
rived directly, clarifying the precise form of time-reversal 
symmetry in the continuum theory. 

The first stage is to express the partition function for 
the bosonic problem as a coherent-state path integral 
(with the hard-core constraint enforced by an explicit 
term added to the Hamiltonian). The quadratic hop- 
ping term is then decoupled using a field ■0^ defined on 
the sites j of the lattice, as in the standard derivation 
of a continuum theory for the bosonic Hubbard model.!^ 
One can then in principle trace over the bosonic degrees 
of freedom to give an effective action for tp. Since the 
field V' couples linearly to the bosonic field b, the effec- 
tive action is symmetric under phase rotation ip — > tpe^^ 
(for arbitrary real 9) and the condensation transition cor- 
responds to the spontaneous breaking of this symmetry 
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of ij}. 

In the case of the bosonic Hubbard model, where there 
is only a quadratic hopping term, decoupling this term 
leaves a single-site problem, for which tracing over the 
bosonic degrees of freedom is tractable. This is no longer 
the case with correlated hopping terms, but the effective 
action for -0 can still be written in the form 



Vbe 



-s[b,4.] 



(A.l) 



Returning to the Hamiltonian representation, this can be 
rewritten as 



g-Seff[V-] = Xr \ T^exp 



(A.2) 



where 'H{ip) is the Hamiltonian TL in which the quadratic, 
hermitian part of the hopping has been replaced by 
The presence of the time-ordering oper- 
ator Tt- (for imaginary time r) makes this form inconve- 
nient to use, and we instead make the time order explicit 
by writing this as 



Tr lim 



(A.3) 



where St = (3/ L and t^ — {^l — |) 5t. 

The symmetry properties of the effective action S^s 
can be determined from this form of the expression. 
First, consider a symmetry such as TZy, defined in Sec- 
tion II B which takes y —y and exchanges the two 
sublattices, but has no effect on the z-direction. This is 
a symmetry of the Hamiltonian 7i, so that we can write 



TZyTlTZy 



Ti, where TZy here denotes the operator rep- 



resenting the transformation (and so TZy 



Uy). For 



H{ip)i we have TZyH{ip)TZy = H{Tiyip), where TZyip indi- 
cates the result of applying TZy to the field ^/;. Note that 
here ip denotes the value of the field at one particular 
instant, say t^. 



Using Eq. (A.3), we can therefore write 



-SaltlTlylp] 



Tr lim I I ( 



(A.4) 



The power-series definition of the exponential and the 
cyclicity of the trace cause the operators to cancel, so 
that Scs[ip] = Scg[TZyip]. 

The same logic applies to particle-hole symmetry at 
half fiUing, since [H,lph] = and so lphH{ip)lph — 

mr)- 

When the Hamiltonian is not particle-hole symmet- 
ric, however, the latter identity is no longer true. If the 
Hamiltonian Ti were hermitian, we would have 'H{ip*) = 
\H{'ip)]\ but instead taking the hermitian conjugate 
changes the direction of the directed hopping, and so we 
have I'yn{ip*)Ty = [H{Ty^)]^ . We write I'y rather than 
Xy , since ip again denotes the value of the field at one par- 
ticular instant, and so only inverts the y-coordinates. 



The analogue of Eq. (A.3) is now 



^-SMi'yi'] = Xr lim ff e-^«[«('^*(-'))]'^i*- , (A.5) 



where the Xy operators can again be cancelled (and the 
power-series definition of the exponential yet again used) , 
leaving 



-5otf[i;^] ^ Xr lim 



(A.6) 



The property of the hermitian conjugate that B'^ = 
{BAy now allows the product to be reversed; after 
changing the labelling of the time steps, we have finally 
Scs[X'y^j] — Scs[ip*], where ip denotes the time-reverse of 
-0. The operator Xy, which effects time reversal as well as 
y-reflection, then obeys the identity ScttlXyip] — Sosl'ip*], 
as assumed in Section HhI 
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